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On a bounded, measurable domain of non-negative current-quark mass, realistic models of QCD's 
gap equation can simultaneously admit two inequivalent dynamical chiral symmetry breaking 
(DCSB) solutions and a solution that is unambiguously connected with the realisation of chiral 
symmetry in the Wigner mode. The Wigner solution and one of the DCSB solutions are desta- 
bilised by a current-quark mass and both disappear when that mass exceeds a critical value. This 
critical value also bounds the domain on which the surviving DCSB solution possesses a chiral 
expansion. This value can therefore be viewed as an upper bound on the domain within which a 
perturbative expansion in the current-quark mass around the chiral limit is uniformly valid for phys- 
ical quantities. For a pseudoscalar meson constituted of equal mass current-quarks, it corresponds 
to a mass m - ~ 0.45 GeV. In our discussion we employ properties of the two DCSB solutions of the 
gap equation that enable a valid definition of (qq) in the presence of a nonzero current-mass. The 
behaviour of this condensate indicates that the essentially dynamical component of chiral symmetry 
breaking decreases with increasing current-quark mass. 

PACS numbers: 12.38.Aw, 12.38.Lg, 11.30.Rd, 24.85. +p 



I. INTRODUCTION 



Dynamical chiral symmetry breaking (DCSB) is the 
creation, via interactions with the gauge field alone, of 
a fermion mass gap: whose magnitude exceeds, per- 
haps by a great amount, the mass-scale in the action 
set by the fermion's bare mass; and which persists 
when that bare-mass-scale vanishes, namely, in the chiral 
limit. It is fundamentally important in strong interac- 
tion physics. For example, DCSB is responsible for the 
generation of large constituent-like masses for dressed- 
quarks in QCD, an outcome that could have been antic- 
ipated from Refs. [H, Q; it is a longstandingprediction of 
Dyson-Schwinger equation (DSE) studies [3| and has re- 
cently been observed in numerical simulations of lattice- 
reg ularised QCD 0,[|. DCSB is also the keystone in the 
realisation of Goldstone's theorem through pseudoscalar 
mesons in QCD [6], and thereby the remarkably small 
value of the ratio of n- and p-meson masses, and the 
weak 7T7T interaction at low energies 0]. 

A large body of efficacious QCD phenomenology is 
built on an appreciation of the importance of DCSB. 
That is evident in studies based on four-fermion inter- 
action models [1, [t| [l(| EL EH an d in DSE applications 
EE Q El- Nevertheless, not all facets of DCSB have 
been elucidated. Herein we describe novel aspects of the 
interplay between explicit and dynamical chiral symme- 
try breaking. 



II. DYNAMICAL CHIRAL SYMMETRY 
BREAKING 

DCSB can be explored via the gap equation; viz., the 
DSE for the dressed-fermion self-energy, which for a given 
quark flavour in QCD is expressed [l6| 
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where represents a Poincare invariant regularisa- 
tion of the integral, with A the regularisation mass- 
scale [1, E3], D^ v (k) is the dressed-gluon propagator, 
r v (q,p) is the dressed-quark-gluon vertex, and m hm is 
the quark's A-dependent bare current-mass. The quark- 
gluon-vertex and quark wave function renormalisation 
constants, Zi.2(£ 2 , A 2 ), depend on the renormalisation 
point, £, the regularisation mass-scale and the gauge pa- 
rameter. 

The solution of the gap equation can be written in the 
following equivalent forms: 



S(p) = 
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(NB. The mass function, M(p 2 ) = B(p 2 , C 2 )/ ' A{p 2 , C 2 ), is 
independent of the renormalisation point.) It is obtained 
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from Eq. ([T]) augmented by the renormalisation condition 



S{p) 1 \ n 2-t2 = H-P + miO- 



(4) 



where m(C) is the renormalised (running) mass: 

Z 2 (C 2 ,A 2 )m bm (A) = Z 4 (C 2 ,A 2 )m(C), (5) 

with Zi the Lagrangian-mass renormalisation constant. 
In QCD the chiral limit is strictly and unambiguously 
defined by %\r% 



Z 2 (C 2 ,A 2 )m bm (A) =0,VA»C 



(6) 



which states that the renormalisation-point-invariant 
current-quark mass rh = 0. 

QCD's action is chirally invariant in the chiral limit. 
Consider a global chiral transformation applied to one 
particular flavour of quark, characterised by an angle 9. 
Under this operation the quark's propagator is modified: 



S(P) 
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(7) 

Suppose DCSB takes place so that B(p 2 ,C i 2 ) ^ 0. 
Then, with the choice 6 = n/2, Eq. (J7J) corresponds 
to mapping B{p 2 ,Q 2 ) — > —B(p 2 ,( 2 )- It follows that if 
B(p 2 ,( 2 ) is a solution of the gap equation in the chi- 
ral limit, then so is [— B(p 2 , ( 2 )]. While these two so- 
lutions are distinct, the chiral symmetry entails that 
each yields the same pressure [18|. Hence they corre- 
spond to equivalent vacua. This is an analogue of the 
chiral-limit equivalence between the (<r = l,7r = 0) 
and (a = — 1, 7T = 0) vacua in the linear-sigma- model, as 
elucidated in Refs. [HI HTJ. It is notable that, more 
generally, given a solution of the m(£) > gap equa- 
tion characterised by {A m{() (p 2 , C 2 ), B m{Q (p 2 , ( 2 )}, then 

{A- m ((){p 2 , ( 2 ), --B-m(C) (P 2 , C 2 )} is a solution of the gap 
equation obtained with [— m(£)]. 

Studies of DCSB have hitherto focused on a positive 
definite solution of the gap equation because the intro- 
duction of a positive current-quark bare-mass favours this 
solution; viz., if another solution exists, then it has a 
lower pressure. Returning again to the sigma-model anal- 
ogy, such a bare-mass tilts the so-called wine-bottle po- 
tential, producing a global minimum at (er = 1, it = 0). 
However, whether the massive gap equation admits solu- 
tions other than that which is positive definite, the effect 
of the current-quark mass on such solutions, and their 
interpretation, are questions little considered. 



III. EXEMPLAR 

To begin addressing these questions, we first consider 
the simple example defined by Eqs. ([I]), ([2]) with the fol- 
lowing forms for the dressed-gluon propagator and quark- 
gluon vertex: 

g 2 D^(p-q) = S MV -^6(A 3 -(?), (8) 
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(9) 



wherein mq is some "gluon" mass-scale and A serves as 
a cutoff [2l|. The model thus obtained is not renormal- 
isable so that the regularisation scale A, upon which all 
calculated quantities depend, plays a dynamical role and 
the renormalisation constants can be set to one. In the 
model thus defined the gap equation is 
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This gap equation's solution is A(p 2 ) 
B{p 2 ) = M, a constant which satisfies 



(10) 
1 and 



M = 771 
f2 a2\ A 2 ji/f2 



M-L-L^M 2 ,! 2 ), (11) 
3tt z tuq 



C(M 2 ,A 2 ) = A 2 -M z \n 1 + A /M 



(12) 



Since A defines the mass-scale in a nonrenormalisable 
model, we can set A = 1 and hereafter interpret all other 
mass-scales as being expressed in units of A, whereupon 
the gap equation becomes 

G(M) := M - 77i bm - M —5- C*(M 2 , 1) = . (13) 

37r^ m G 



, bm 







Equation (TT3")) admits a M ^ solution when 
if and only if 

m 2 G < (mg) 2 = ^ ; 

namely, it supports DCSB in this case. Hence, to proceed 
we choose 



(14) 
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NB. For me > m c e the only solution of the gap equation 
is one that may be obtained via a perturbative expansion 
in the coupling and hence DCSB is impossible. That 
domain is therefore not of interest herein. 

G(M) is plotted in Fig.Q] One reads from the figure 
that in the chiral limit there are three solutions to the 
gap equation: 



M = 



M w = , 

M± = ±M° = ±0.33 . 



(16) 



M\y describes a realisation of chiral symmetry in the 
Wigner mode. It corresponds to the vacuum configura- 
tion in which the possibility of DCSB is not realised; i.e., 
a = = 7T in the sigma-model analogy. In the chiral 
limit this is the only solution accessible via a perturba- 
tive expansion in the coupling. The solutions M± are 
essentially nonperturbative. They represent the realisa- 
tion of chiral symmetry in the Nambu-Goldstone mode; 
namely, DCSB. 




FIG. 1: The zeros of G(M) give the solution of the gap equa- 
tion defined by Eqs. (|13|l . (|15[) . Solid curve: obtained with 
m bm = 0, in which case G(M) is odd under M -> —A/; 



long-dashed curve: m 
m c ™ = 0.033; dotted curve: 
quantities in units of A.) 



0.01; short-dashed curve: m = 
' m = 0.05. (All dimensioned 



It is apparent in Fig. [T] that with the gap equation ob- 
tained via Eqs. ©, ©, each of the solutions identified in 
Eq. p6p evolves smoothly with current-quark mass on a 
neighbourhood of m bm = 0. We subsequently consider 
the manner in which these solutions evolve as m bm is in- 
creased from zero. In that discussion we will retain the 
labels introduced in Eq. (fl6|) and attach them to that 
solution which is (pointwise, if relevant) closest in mag- 
nitude to the chiral limit solution of the same name. 

M = M + > is the solution usually tracked in connec- 
tion with QCD phenomenology. In models of this type 
it is identified as a constituent-quark mass. As m bm is 
increased, M + also increases. 

As evident in Fig.[Tl the other two solutions of the 
gap equation do not immediately disappear when m hm 
increases from zero. Nor do they always persist. Instead, 
these solutions exist on a domain 



V(m hm ) = {r 



< m bm < m b r m }. 



(17) 



JV%,M_ also evolve smoothly with w bm . Moreover, at 
the critical current-quark mass, m b ™, these two solutions 
coalesce. 

To understand the origin of a critical mass we observe 
from Eq. (|13[) and Fig.[T]that introducing a current-quark 
mass merely produces a constant pointwise negative shift 
in the curve G{M). Hence, the critical current-quark 
mass is that value of this mass for which G(M) = at 
its local maximum. The local maximum occurs at 



M lm = {M \M < , G'(M)} = 



(18) 



and therefore in the present illustration the critical 
current-quark mass 

m bm = {m hm | G(M\ m ) = 0} = 0.033 . (19) 



FIG. 2: Circles - Reciprocal of the n-th roots of the coeffi- 
cients in a chiral expansion of M+ around m bm = 0, Eq. (|20|l . 
as a function of t — l/n. Dashed curve - the function in 
Eq. |H3]|. Dotted curve - m b r m = 0.033. (All dimensioned 
quantities in units of A.) 



We now return to the behaviour of M + and ask 
whether this quantity has a power series expansion in 



in 



bin 



about m 



bm 



0; viz., a chiral expansion: 



M+(m bm ) =Mo + Y a n (m hnl Y 



(20) 



Such an expansion exists; i.e., is absolutely convergent 
on a measurable domain, so long as V n 



1 



l/n 



> 



1 



!/("+!) 



and 



l/n 

:= lim I — — I > , 



(21) 



(22) 



where the quantity m rc is the radius of convergence for 
the series; i.e., the series converges on < m bm < m rc . 

To determine whether M + has such an expansion we 
inserted Eq. (|20|) into Eq. (fT"3"f and solved the sequence 
of algebraic equations that this produces to obtain the 
coefficients {a\, . . . ,a 14 }. The procedure is straightfor- 
ward but we stopped at n — 14 because the mag- 
nitude of the coefficients grows rapidly with n; e.g., 
ai4 = —4.27331 x 10 17 , and this order was sufficient for 
our purpose. The information depicted in Fig.[2]indicates 
that M_|_ does have a chiral expansion with a nonzero but 
finite radius of convergence. The curve in the figure is 
the function (m = 1 diagonal Pade) 



U\t 



Uq 



1 + u 2 t 



, it = 0.326 , ui = 0.295 , u 2 = -0.596 , 

where the coefficients were fixed in a least-squares fit to 
{(l/a n ) 1/n , n = 2, . . . , 14}. The limit t -> corresponds 
to n — > oo and hence the value of uq gives the radius 



4 



0.3 



0.2 



0.1 



_i i i_ 



_i i i i_ 



_i i i i_ 



0.01 



m 



0.02 

bm 



0.03 



FIG. 3: Evolution with current-quark mass of M (solid curve) 
and M (dashed curve). (All dimensioned quantities in units 
of A.) 



of convergence. We repeated the fit with diagonal Pade 
approximants of order m — 2,3,4. Combining the re- 
sults, we obtain the radius of convergence for the chiral 
expansion of M + : 



0.034 ±0.001 



.bm 



(24) 



within numerical error. Hence D(m bm ) is also the domain 
on which M + has a chiral expansion. (NB. Analogous 
chiral expansions exist within this domain for Mw, M-. 
The expansion for Mw vanishes for m bm — whereas 
the leading term in that for M_ is Mq. We emphasise 
that M± have chiral expansions around the solution that 
is nonperturbative in the coupling. A nonzero value of 
Mq is never attainable via an expansion in the coupling.) 
The two combinations 



M 



1 



(M 4 



M_) , M 



1 



(M++M_) (25) 



are of interest. In the chiral limit, M = M° and 
M = Mw- The evolution of each with current-quark 
mass is depicted in Fig. [31 It is apparent that M is con- 
tinuous on X>(ra bm ) and evolves from the DCSB solu- 
tion with increasing m : M(m bm ) is a monotonically 
decreasing function. These features can be understood 
as illustrating that the essentially dynamical component 
of chiral symmetry breaking decreases with increasing 
current-quark mass. This has also been argued via the 
constituent-quark cr-term; e.g., in Sec. 5.2.2 of Ref. (22J. 

The alternative combination, M, is also continuous 
on T>(m hm ). With the value of the coupling given in 
Eq. (|15p. M(m bm ) evolves from the Wigner solution ac- 
cording to 



M m = 



2 3(M°) 2 - 1 
1+ 3 (M°) 2 + l 



(26) 



The development of M(m bm ) might be viewed as a gauge 
of the destabilising effect that DCSB has on this model's 
Wigner mode. 

We now return to the critical current-quark mass, 
Eq. fig]). (j2"4"]l . In the neighbourhood of m bm = 0, 



M w (m hm ) = -3m bm + 



(27) 



viz., Mw(m ) can be expressed as a power series in 
m hm around its chiral limit value, where that value is 
perturbative in the coupling. However, with increasing 
current-quark mass Mw decreases steadily toward M_, 
which is nonzero in the chiral limit and essentially non- 
perturbative in the coupling, until at m b ™, Mw — M_. 
At this point a solution whose small current-quark mass 
behaviour is essentially perturbative has melded with a 
solution that is inaccessible in perturbation theory and 
actually characteristic of DCSB. A related view sees M_ 
as a DCSB solution whose modification by a current- 
quark mass can no longer be evaluated as a power series 



in m when that mass exceeds 



Finally, a chi- 



ral expansion of M+ can only converge for current-quark 
masses less than m b ™. 

These observations suggest that m b ™ specifies the up- 
per bound on the domain within which, for physically rel- 
evant quantities, a perturbative expansion in the current- 
quark mass around their chiral limit values can be valid; 
i.e., it is a (possibly weak) upper bound on the radius 
of convergence. This view and the numerical result in 
Eq. (fTO)) coincide with Ref. [23| . We will return to this 
point. 



IV. CLOSER TO QCD 

It is natural to ask whether analogous behaviour exists 
in QCD. To explore this we work with a renormalisation- 
group-improved (RGI) rainbow truncation of the gap 
equation's kernel. This is the leading-order in a system- 
atic and symmetry-preserving truncation of the DSEs 
that is noperturbative in the coupling [24], HH, 1261 l27l . 
The truncation has been used widely; e.g., Refs. [17], [H, 
[29j . and references thereto, and an efficacious implemen- 
tation preserves the one-loop ultraviolet behaviour of per- 
turbative QCD. However, a model assumption is required 
for the behaviour of the kernel in the infrared; viz., on the 
interval Q 2 < 1 GeV 2 , which corresponds to length-scales 
> 0.2 fin. 

The rainbow truncation is realised in the gap equation 
via the replacement of Eq. @ by (TtJ 



A" 



m = J G{(V-<lYW*7(p-<l)^l»S{q)- lv . (28) 
Herein we employ the model interaction introduced in 
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Ref. [30|1: 



G(t) _ 4tt s 



Dte 



-t/u 



8tt 2 



In 



t/A 2 



QCD 



(29) 



(4) 
'MS 



with £ = fc 2 , J"(t) = [1 - cxp(-t/[4m 2 r ])]/i, m 2 , 
0.5 GeV, t = e 2 - 1, 7m = 12/25 and A QCD 
0.234. This form expresses the interaction as a sum of two 
terms. The second ensures that perturbative behaviour 
is correctly realised at short range; namely, as written, 
for (k - q) 2 ~ k 2 ~ q 2 > 1 - 2 GeV 2 , Eq. (22) guar- 
antees that the quark-antiquark scattering kernel, AT, is 
precisely as prescribed by QCD. On the other hand, the 
first term in Q (t) is a model for the long-range behaviour 
of the interaction. It is a finite width representation of 
the form introduced in Ref. [31( , which has been rendered 
as an integrable regularisation of 1/t 2 = 1/k 4 [32| . This 
interpretation, when combined with the result that in a 
heavy-quark-heavy-antiquark BSE the RGI ladder trun- 
cation is exact [23|, is consistent with Q{t) leading to a 
Richardson-like potential [H| between static sources. 

The true parameters in Eq. (|29|) are D and uj, which 
together determine the integrated infrared strength of the 
rainbow kernel. However, they are not independent [30j | : 
in fitting to a selection of ground-state observables, a 
change in one is compensated by altering the other; e.g., 
on the interval uj € [0.3, 0.5] GeV, the fitted observables 
are approximately constant along the trajectory [34] 



luD = (0.72 GeV) 3 



(30) 



Equation (|29p is thus a one-parameter model. 

It is important to bear in mind that because the 
truncation preserves the one-loop renormalisation group 
properties of QCD the ultraviolet behaviour of the solu- 
tions of Eqs. ([1]) and (j2"5|) is precisely that of QCD. Hence 
we have 



Af(p 



2 x P »Aqcd 



(IlnpVA 2 QCD 



(31) 



where rh is the renormalisation-group-invariant mass. 

The model-dependence is mainly restricted to infrared 
momenta but on this domain, too, there is good agree- 
ment with QCD; e.g., the gap equation solutions are 

with numerical sim- 



34] 



in semiquantitative agreement 
ulations of lattice-regularised quenched-QCD. (NB. Pre- 
cise agreement would be incorrect because Eq. (|2"9"|) cor- 
responds to an unquenched theory.) The conditions have 
been explored under which pointwise agreement between 
DSE results and quenched- and full-lattice simulations 
may be obtained [Hi H H3 . 

Equations ^ and (|2"5j) obviously admit the M(p 2 ) = 
solution in the chiral limit, Eq. ©. This solution, which 
can always be obtained through a weak coupling expan- 
sion, is analogous to Mw in Eq. P^|) . 

In the chiral limit the rainbow gap equation also yields 
a DCSB solution. This capacity is the basis for much of 




0.01 0.02 0.03 

m (c 19 ) [GeV] 



0.04 



FIG. 4: Evolution with current-quark mass, m((ig), of 
A(p 2 = 0,( 2 ) (dimensionless), B(p 2 = 0,C 2 ) (GeV) as cal- 
culated with uj = 0.4 GeV in Eq. ([29|l : solid curve - B+(0); 
dash-dot curve - A+(0); long dashed curve - B-(0); sparse 
dotted curve - A-(0); short dashed curve - Bw(0); dense 
dotted curve - Aw(0)- 



the phenomenological success of the RGI rainbow-ladder 
truncation. The truncation preserves the feature that if 
M + (p 2 ) = M(p 2 ) > 0, Vp 2 > 0, is a solution of the chiral- 
limit gap equation, then so is M_(p 2 ) := [— M(p 2 )]. 
These solutions are the analogues of M± in Eq. (fT6f . 

We solve the gap equation using the Pauli-Villars reg- 
ularisation procedure described in Ref. 38] and work 
in Landau gauge because it is a fixed point of the 
renormalisation group. The renormalisation condition, 
Eq. (J4|, is implemented at £ = fi9 := 19 GeV with a 
choice for the current-quark mass. This scale is deep 
in the perturbative (weak-coupling) domain and hence 
a renormalisation-group-invariant current-quark mass is 
unambiguously specified via Eq. (|31[) . For reference, we 
also list below current-quark masses defined therewith 
via one- loop evolution to a "typical hadron scale"; viz., 



m(Ci 



(lnCi/AqcD) 7 " 



Ci = 1 GeV . 



(32) 



It is apparent in Fig. 3] that this interaction model, 
too, exhibits a bounded domain of current-quark mass 
on which M_(p 2 ) and M + (j> 2 ) exist simultaneously: 



T>(m) — {rh \ < m < m cr } 



(33) 



Away from the chiral limit the solution characterised 
by M_(p 2 ) describes the propagation characteristics of 
a quark embedded in an unstable vacuum. The prop- 
erties of that vacuum may be detailed, e.g., by solv- 
ing the Bethe-Salpeter equation for meson bound-states 
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i 1 — i — i i i 1 1 1 1 1 — i — i i i 1 1 1 1 1 — i — i i i 1 1 





10~ 2 10~ 1 10° 10 1 
p 2 [GeV 2 ] 



FIG. 5: Momentum dependence of the dressed-quark mass- 
function, M(p 2 ): upper panel - chiral limit; middle panel - 
m(Cig) = 5MeV; lower panel - m(£i9) = 50MeV, for which 
there is naturally no M-(p 2 ) solution. In each panel the solid 
curve is M+(p 2 ); the dashed curve is M_(p 2 ); and the dotted 
curve is Mw(p 2 )- All results obtained with uj = 0.4 GeV. 



with this dressed-quark propagator. The critical current- 
quark mass depends on uj: 



uj [GeV] 


0.3 


0.4 


0.5 


m cr [MeV] 


71 


63 


31 


m cr (Ci 9 ) [MeV] 


35 


31 


15 


TOcr(Cl) [MeV] 


60 


53 


26 



(34) 



where the third and fourth rows report the critical mass 
at the renormalisation scales described above. (NB. In 
most phenomenological applications 0.3 < uj < 0.4.) In 
analyses of hadron observables founded on the models 
that form the basis of our arguments, the s-quark RGI 
current-mass exceeds m cr . 

It is noteworthy that we also find a Wigner solution 
when applying Eq. ([J); i.e., in this case, for nonzero 
current-quark mass we find an analogue of Myv in 
Eq. (|16p . This suggests that the presence of the Wigner 



solution is not contingent upon the pointwise behaviour 
of the gap equation's kernel and supports an interpreta- 
tion of our findings in the context of QCD. (Indeed, all 
of the features identified herein are also expressed, e.g., 
in the model of Ref. [39(.) We observe that in the model 
of this section the Wigner solution is naturally repre- 
sented by two momentum-dependent functions; namely, 
A w (p 2 X 2 ), B w (p 2 X 2 )- As indicated by Figs.H El E 
at rh = rhcx these functions meld pointwise with their 
counterparts in the DCSB M_(p 2 ) solution. 

We emphasise that all solutions of the gap equation 
evolve smoothly with current-quark mass on the domain 
specified in Eq. (f34j) . Moreover, the Nambu solution char- 
acterised by M + (p 2 ) evolves smoothly for all values of 
the current-quark mass. Naturally, this does not nec- 
essarily entail that it's pointwise behaviour at a given 
value of current-quark mass is obtainable via a perturba- 
tive expansion in rh about the nonzero DCSB result at 
to = 0, which is essentially nonperturbative in the cou- 
pling. Indeed, the commonality of behaviour between the 
two models we discuss explicitly herein, and with the oth- 
ers we have considered, suggests strongly that m cr is the 
radius of convergence for a pointwise chiral expansion of 
M + (p 2 ) around the DCSB chiral-limit result. However, 
verification must wait because the method employed in 
arriving at Eq. (|2"4]) is not workable for the integral gap 
equation and we do not yet have a tractable alternative. 

In considering other rainbow interaction models, we 
find that to ci . assumes similar values in most models that 
provide a reasonable description of the same low-energy 
observables. An exception is the model of Ref. [3l|, in 
which a solution with B(s — 0) < 0, and A(s) and B(s) 
continuous on s G [0, oo) exists only in the chiral limit, 
with the interpretation that in this model to„. = 0. 

One might also ask after the effect of dressing the 
quark-gluon vertex. As noted above, symmetry ensures 
that in the chiral limit the gap equation simultaneously 
admits M_(p 2 ) and M + (p 2 ) solutions in this instance, 
too. The extent of the domain of current-quark mass on 
which the M_(p 2 ) solution persists will probably depend 
on the structure of the vertex. We are currently exploring 
this. 

Some may seek an interpretation of the gap equation's 
distinct solutions. Here, again, the tr-model analogy is 
useful. Consideration reveals that the gap equation's so- 
lutions are representations at an elementary level of the 
misaligned chiral order parameter that is the focus in 
discussions of disordered chiral condensates, many con- 
sequences of which are reviewed in Ref. [40( • We are cur- 
rently examining others that may be equally useful; e.g., 
as stated above, by calculating hadron properties in the 
different vacua. 

Equation is valid for M±(p 2 ) and M w (p 2 ). In 

fact, one can make a stronger statement, while M±(p 2 ), 
M w {p 2 ) exist: 

M+ (p 2 ) p2>> ^ CD AUp 2 ) p2>> = QCD M W ( P 2 ) . (35) 
This follows from asymptotic freedom, which is a feature 
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FIG. 6: Momentum dependence of the dressed-quark mass- 
function: dashed curve - M_(p 2 ); dotted curve - Mw(p 2 ), 
both obtained with m(( 19 ) = 5 MeV and u = 0.4 GeV. 



FIG. 7: Evolution with current-quark mass, m((ig), of the 
massive-quark condensate defined in Eq. (|37f) . calculated with 
uj = 0.4 GeV and at the renormalisation point £ = 19 GeV. 



of the RGI model and QCD. One may argue for this result 
as follows. On the weak-coupling domain 



A±w(p 2 X 2 ) « 1 



M± tW (p 2 



P 2 + MLw(P 2 ) 



M±,w(p 2 



(36) 



Hence, the gap equation becomes a single linear integral 
equation for M± w(p )■ Within the domain on which 
the preceding steps are valid, that equation can be ap- 
proximated by a linear second-order ordinary differential 
equation (d.e.) (e.g., Refs. [M S3, SI S3)- The d.e. is 
the same for every one of the functions M± j w{p 2 ), as is 
the ultraviolet boundary condition, which is determined 
by the current-quark mass. Thus follows Eq. (|35p . a re- 
sult evident in Figs. 03 [H 

In general the gap equation is solved as a nonlinear 
integral equation. It is straightforward to obtain the 
M±(p 2 ) solutions via iteration. However, for nonzero 
current-quark mass the Wigner solution is harder to fix. 
A careful seed function for iteration must be chosen and 
an adaptive iterative approach employed to reach the so- 
lution. To be concrete, a Chebyshev expansion was em- 
ployed for the solution functions and a Newton iteration 
procedure used to determine the coefficients, after the 
manner of Ref . [45] . The seed for iteration was A = 1 
and B = — 10m(£), and convergence to solution first ob- 
tained with loD = (0.6 GeV) 3 . With uj = 0.4 GeV, small 
steps in D were subsequently made to reach the value in 
Eq. (|30[) . All solutions are illustrated in Fig. [5] and, for 
clarity and emphasis, we compare M_(p 2 ) and Mw{p 2 ) 
in Fig.[6l 

In the context of the d.e. argument presented above, we 
add that the model gap equation has often been approx- 
imated by a single second-order d.e., which is nonlinear 



in M(p 2 ) for infrared momenta but linear in the ultra- 
violet [H SI SI, S3- As 

we remarked, the ultraviolet 
boundary condition for all solutions is still fixed by the 
current-quark mass and the solutions agree. However, as 
apparent in Fig.[U A\y ^ 1 on the domain of infrared mo- 
menta. Hence, while a single d.e. remains a valid approx- 
imation for M±(p 2 ), that is not the case for the Wigner 
solution. This emphasises that the differences between 
M + (p 2 ), M_(p 2 ) and Mw(p 2 ) are a primarily infrared 
effect; i.e., nonperturbative in the coupling. 

Using Eq. (|55|) and the subsequent discussion it be- 
comes clear that the gauge-invariant quantity 



<7(m(C)):= lim Z A (( 2 , A 2 ) N c tr D / 5 m (%,0, (37) 



where 



S m{C) (q,C) = g 



(38) 



is a current-quark-mass-dependent quark condensate 
that is well-defined, finite and unambiguous, and has 
a perturbative expansion in rh on a bounded interval. 
In addition, tr(m(C)) evolves under the renormalisation 
group in precisely the same manner as the chiral-limit 
vacuum quark condensate and is identical to the vacuum 
quark condensate in the chiral limit. Equation (|37| is 
unique in possessing all these properties. 

For additional clarity we remark that for rh ^ there 
is no term in a perturbative expansion around the chiral 
limit of the integrand in Eq. ([3T)l that is odd in rh. Now, 
since a weak-coupling evaluation of the scalar piece of 
the quark propagator yields an expression that is odd 
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in TTi, then Eq. (|37|) contains no term calculable in weak- 
coupling perturbation theory. (This may be verified using 
the d.e. analysis discussed above.) 

The behaviour of a(m(Q) obtained with the RGI 
model employed herein is depicted in Fig. [71 As we 
saw in connection with Fig. [3J here, too, the essen- 
tially dynamical component of chiral symmetry break- 
ing decreases with increasing current-quark mass, follow- 
ing the trend predicted by the constituent-quark cr-term 
(Ref. [H, Sec. 5.2.2). 

At this point we can further amplify our interpreta- 
tion of m cr . A basic ingredient of chiral perturbation 
theory is the quark condensate because it introduces the 
mass-scale of DCSB. If the evolution with current-quark 
mass of this quantity about its DCSB chiral-limit value 
cannot be evaluated perturbatively then the expansion 
in to around to = has broken down. The pseudoscalar- 
meson-loop contribution to this evolution is a long-range 
piece of the vacuum response. Now, with Eq. (|37[) . we 
can evaluate a short-range piece. While an expansion of 
both M- and M+ in to exists, then so does a pertur- 
bative expansion of the order parameter a. However, a 
has no such expansion for to > m cr . NB. A chiral ex- 
pansion is invalid if either the long-range or short-range 
contribution fails but it is likely that they are linked. 

We stress that a straightforward definition of a 
massive-quark condensate via the trace of a to 7^ 
dressed-quark propagator is not useful because it gives 
a quantity that is quadratically divergent and therefore 
very difficult to define unambiguously. The same weak- 
ness afflicts the quantity 

,A 

<t(to(C),A) = Z 4 (C 2 ,A 2 )AT c tr D / S m <%,C), (39) 

Jq 



where 



s m( %, C) = ^ C) + s m( %, C) 



(40) 



The discussion in Refs. p,|46|,|47| provides a context for 
Eqs. (|37| - (|40| . and a connection between our reasoning 
and that used in other frameworks to calculate the quark 
condensate. 



V. SUMMARY 

On a bounded interval of current-quark mass, T>(m) — 
{to I < to < TO cr }, realistic models of QCD's gap 
equation can simultaneously admit two inequivalent dy- 
namical chiral symmetry breaking (DCSB) solutions for 
the dressed-quark mass- function, M±(p 2 ), and a solu- 
tion that is unambiguously connected with the realisa- 
tion of chiral symmetry in the Wigner mode, M\y{p 2 )- 
The DCSB solutions are distinguished by their value at 
the origin: M + (p 2 = 0) > and Af_(p 2 = 0) < 0. In the 
ultraviolet all three solutions coincide with the running 
current-quark mass. 



The pointwise values of all solutions evolve continu- 
ously with current-quark mass within T>(rh). However, 
things change at the upper boundary. The Mw solution, 
whose chiral limit value is perturbative in the coupling, 
becomes identical to the essentially nonperturbative so- 
lution M_ that is actually characteristic of DCSB. More- 
over, both disappear for m > m cr , a domain whereupon 
the current-quark mass is large enough to completely 
destabilise these solutions. Only the positive M + solu- 
tion exists on this domain. Furthermore, we provided ev- 
idence that the upper boundary of X>(to) also defines the 
radius of convergence for an expansion of M + in current- 
quark mass around its DCSB chiral-limit form. 

Thus one has the coalescing of two qualitatively dis- 
tinct solutions at rh CI , the persistence of only one es- 
sentially nonperturbative solution for to > rh CI plus the 
breakdown of a chiral expansion for this solution, and 
the simultaneous loss of a well-defined mass-dependent 
quark condensate. This behaviour supports an interpre- 
tation of the upper bound on the domain within 
which a perturbative expansion of physical quantities in 
the current-quark mass around their chiral-limit values 
can uniformly be valid. In a phenomenologically effi- 
cacious renormalisation-group-improved rainbow-ladder 
truncation of QCD's Dyson- Schwinger equations, the 
critical current-quark mass corresponds to a mass m - ~ 
0.45 GeV (to 2 _ ^0.20GeV 2 ) for a pseudoscalar meson 
constituted of equal mass current-quarks 48]. NB. Irre- 
spective of the current-mass of the other constituent, a 
meson containing one current-quark whose mass exceeds 
TO cr is never within the domain of uniform convergence. 

A value of similar magnitude was deduced in Refs. [49|, 
[50j as the scale below which accuracy may be expected 
from the approximation of observables through a pertur- 
bative expansion in pion-like pseudoscalar-meson mass. 
This scale also marks the boundary below which observ- 
ables should exhibit that curvature as a function of pion- 
like pseudoscalar-meson mass which is characteristic of 
chiral effective theories. 
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